Automatic Segmentation of Pulmonary Artery (PA) Using Customized Level Set Method in 3D (CTA) Images

This paper proposes an efficient customized level set algorithm to segment Pulmonary Artery (PA) tree in 3D pulmonary Computed Tomography Angiography (CTA) images. In this algorithm, to reduce the search area the lung regions from the original image are first segmented and the heart region is extracted by selecting the regions between the lungs. A pre-processing algorithm based on Hessian matrix and its eigenvalues is used to remove the connectivity between the pulmonary artery and other nearby pulmonary organs. To extract the pulmonary artery tree, we first use a region growing method initialized by a seed point which is automatically selected within the pulmonary artery trunk in the heart region. In the second step, the segmentation of the pulmonary artery is performed using a customized 3D level set algorithm, using the output of region grower as the initial contour. We use a new stopping criterion for the proposed level set algorithm, a consideration often neglected in many level set implementations. To validate and assess the robustness of the method, 20 CT angiography datasets were used (10 free pulmonary embolism scans and 10 CT with pulmonary emboli). A very good agreement with the visual judgment was obtained in both normal and positive pulmonary emboli CT scans.